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We study one-dimensional Su-Schrieffer-Heeger (SSH) models with quantum phonons using a 
continuous-time quantum Monte Carlo method. Within statistical errors, we obtain identical results 
for the SSH model with acoustic phonons, and a related model with a coupling to an optical bond 
phonon mode. Based on this agreement, we first study the Peierls metal-insulator transition of the 
spinless SSH model, and relate it to the Kosterlitz-Thouless transition of a spinless Luttinger liquid. 
In the Peierls phase, the spectral functions reveal the single-particle and charge gap, and a central 
peak related to long-range order. For the spinful SSH model, which has a dimerized ground state 
for any nonzero coupling, we reveal a symmetry-related degeneracy of spin and charge excitations, 
and the expected spin and charge gaps as well as a central peak. Finally, we study the SSH-f/H 
model with electron-phonon and electron-electron interaction. We observe a Mott phase with critical 
spin and bond correlations at weak electron-phonon coupling, and a Peierls phase with gapped spin 
excitations at strong coupling. We relate our findings to the extended Hubbard model, and discuss 
the physical origin of the agreement between optical and acoustic phonons. 

PACS numbers: 71.10.Hf, 71.10.Pm, 71.30.+h 


I. INTRODUCTION 

In 1979, Su, Schrieffer and Heeger (SSH) introduced a 
simple microscopic model of polyacetylene with a modu¬ 
lation of the electronic hopping as a result of distortions 
of the carbon bonds [1]. It has since gained a much wider 
relevance due to the fact that it supports topological ex¬ 
citations (solitons) with fractional charge [1-3], and its 
relation to one-dimensional topological insulators [4] in 
the BDI class [5]. As a model for correlated electrons, 
the case of a half-filled band has attracted particular in¬ 
terest, because it supports a Peierls and (in the presence 
of Coulomb interaction) a Mott insulating phase. If the 
lattice is treated quantum mechanically, the model con¬ 
stitutes a rich and complex many-body problem. 

In the adiabatic limit (or at the mean-field level), a 
dimerized Peierls state results at T = 0 for any nonzero 
electron-phonon coupling [1, 6]. Although a classical 
treatment of the lattice (see, e.g.. Refs. [7-9]) appears 
justified for polyacetylene given the large oscillator mass 
[1], quantum fluctuations of the lattice cause a signifi¬ 
cant reduction of the dimerization [6, 10-13] and affect 
the low-lying excitations [14]. Quite surprisingly, numer¬ 
ical [6, 10, 15] and analytical [6, 16-18] results for the 
SSH model without electron-electron interaction suggest 
that lattice fluctuations do not destroy the Peierls state, 
in contrast to the spinless SSH model [6] or Holstein mod¬ 
els [19-21]. Interestingly, the Peierls state does become 
unstable with respect to quantum fluctuations if addi¬ 
tional Coulomb repulsion between electrons is taken into 
account, as appropriate for polyacetylene [22]. In fact, 
the single-particle gap found experimentally has been ar¬ 
gued to be predominantly a result of electronic correla¬ 
tions, with only a small component from the dimerization 
of the lattice [7]. A finite phonon frequency may also 
crucially affect the balance between interactions [23, 24]. 
The SSH models are also related to spin-Peierls models 


[6, 25], studied in detail to understand the physics of 
CuGeOg [26], see for example Refs. [27, 28]. For a more 
complete overview of previous work on polyacetylene we 
refer to Refs. [22, 29]. 

The momentum dependence of the phonon dispersion 
and the interaction make numerical studies of the SSH 
model more challenging than for models with a Holstein 
interaction [30]. Therefore, several open problems re¬ 
main, in particular the calculation of momentum-resolved 
excitation spectra. Besides, some of the recent numeri¬ 
cal work has focused on a simplified model with optical 
phonons and a different interaction term [25]. While the¬ 
ory suggests that the low-energy physics depends solely 
on the phonon frequency at g = 2k^ [6, 16], the quantita¬ 
tive impact of the phonon dispersion is not fully under¬ 
stood. In particular, recent functional renormalization 
group results for the model with acoustic phonons [18] 
agree quantitatively with numerical results for the model 
with optical phonons [25]. On the other hand, a signifi¬ 
cant effect of the dispersion was reported in Ref. [31]. 

Here, we use quantum Monte Carlo simulations of the 
original SSH model (with acoustic phonons) to demon¬ 
strate that it gives results that are identical within statis¬ 
tical errors to the simplified model with optical phonons 
for all parameters considered. Given this agreement, we 
study in detail the real-space correlation functions and 
calculate spectral functions of the simplified models. We 
consider the spinless and the spinful SSH model, as well 
as the SSH model with additional Coulomb interaction. 

The paper is organized as follows. In Sec. II we define 
the models, in Sec. Ill we discuss the key methodological 
aspects, in Sec. IV we present the numerical results, in 
Sec. V we discuss the impact of the phonon dispersion 
and lattice fluctuations as well as the low-energy theory, 
and in Sec. VI we conclude. 
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II. MODELS 


C. SSU-UV model 


In addition to the original SSH model [1], two other 
variants have been studied in the past. First, an SSH 
model with optical phonons but with the original cou¬ 
pling term, proposed in Ref. [32], and second a model 
with optical phonons and a simplified coupling term in¬ 
troduced in Ref. [25]. 


A. SSH model 


Defining the bond operator 

^ V ^^icr^i + lcr ^i+lcr^icr^ (^) 

a 


the SSH Hamiltonian [1] can be written as 


( 2 ) 


E 


—P^ + — (Q , 1 - Q p 
2M " 2 ^ 


The first term describes nearest-neighbor electronic hop¬ 
ping with amplitude t in the undistorted chain, the sec¬ 
ond term is a coupling between electrons and phonons in 
the form of a modulation of the electronic hopping as a 
result of lattice distortions (see also Ref. [33]), and the 
last term describes harmonic oscillators with mass M, 
spring constant AT, coordinate Q-, and momentum P-; it 
can be written as with the phonon dispersion 

Wg = w^sin(g/2). 

The spinless SSH model has the same form as Eq. (2), 
with the bond operator given by B- = [^]* 


B. Optical SSH model 


We also simulated the Hamiltonian [25] 




P— p2 ^ —Q'^ 
2M * 2 




(3) 


which describes the coupling o f elec trons to optical 
phonons with frequency uJq = ^JK/M. As in the SSH 
model, lattice distortions modulate the electronic hop¬ 
ping, but the interaction only involves the coordinate Q- 
rather than the difference ~ Qi' following, 

we refer to Eq. (3) as the optical SSH model, and will 
demonstrate that it gives the same results as the SSH 
model (2). The optical phonons in Eq. (3) may be re¬ 
garded as describing fluctuations of the bonds rather than 
atom positions as in Eq. (2). As discussed in Sec. V, this 
interpretation provides a physical understanding of the 
observed agreement between Eqs. (2) and (3). 


In the spinful case, we will also study the impact of 
electron-electron interactions. Eollowing previous work 
[7, 23, 25, 31, 34], we considered onsite and nearest- 
neighbor repulsion, as described by 

Ke = + - 1) ’ 

i i 

(4) 

with ft- = We have simulated Eq. (4) in combi¬ 

nation with the original SSH model (2) (SSH-I/E model) 
and the simplified model (3) (optical SSH-I/E model). 


D. Spin-charge symmetry 

Eor our choice of zero chemical potential /i = 0, and on 
bipartite lattices, the models (2) and (3) are invariant un¬ 
der the particle-hole transformation i-^ (— 
and hence half-filled. Eor 1/ = R = 0, the spinful SSH 
and optical SSH models are furthermore invariant under 
[19], which interchanges spin and charge 
degrees of freedom. Whereas this transformation maps 
the repulsive to the attractive Hubbard model [35], here 
it implies a perfect degeneracy of spin and charge degrees 
of freedom for 1/ = E = 0, as confirmed by our results. 


E. Conventions 

We focus on half-filling, corresponding to an average 
density {ft-) = 1 for the spinful case, and to {fi-) = 0.5 
for the spinless case. We use t as the unit of energy, 
and set the lattice constant and h to one. We define the 
dimensionless coupling constant A = /Kt^ and consider 
chains with L sites and periodic boundary conditions. 
The inverse temperature was chosen as = L. 


HI. METHOD 

We used the continuous-time interaction-expansion 
(CT-INT) quantum Monte Carlo method, which is based 
on a weak-coupling expansion of the partition function, 
and gives exact results for finite systems [36]. 

We write the SSH model (2) in the generic form 

^ = -^0 + E + 7 + p\bq) (5) 

q q 

with the charge operator 

Pq = ^ A - E ^kJk+qa [sin k - sm{k + q)] . (6) 
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The partition function takes the form 


A. Sign problem 


Z = y y . (7) 

Here, Sq describes the noninteracting fermions, whereas 
the phonon and electron-phonon contributions to Eq. (2) 
are contained in Carrying out the Gaussian inte¬ 

grals over the phonon degrees of freedom [37, 38] gives a 
retarded and nonlocal electron-electron interaction 


Si = d-rdr y^B^(r)r(z- 

( 8 ) 

where r(r, t) = L ^ e r) and 


D{q,T) 


CQsh[a;g(|T| -^/2)] 

2 smh[uj^l3/2] 


(9) 


A similar representation of the phonon-mediated interac¬ 
tion was given in Ref. [6]; however, in the continuum ap¬ 
proximation, the SSH interaction becomes local in space, 
as characteristic of Holstein models [30]. 

For the optical SSH model, the interaction reduces to 


5i = -^ drdr ^ B^{T)D{T - r )-B.(r'), (10) 


where D{t) follows from Eq. (9) by replacing oOq. 

Although both bond operators in Eq. (10) carry the same 
index i, the interaction remains nonlocal because B.{t) 
involves the sites i and i + 1. 

The interactions described by Eqs. (8) and (10) can be 
simulated with the CT-INT method by implementing a 
general vertex with four lattice indices, two time indices, 
and two spin indices. For the optical SSH models, there 
is no sign problem so that a sampling over additional 
Ising spins [38] is not necessary. 

In the case of the SSH-[/E model, the additional 
Coulomb interaction described by Eq. (4) enters the par¬ 
tition function in the form of 

S2 = U f dr^K^(r) - i][n^^(T) - i] (11) 

Jo ^ 

+V f dr ^ - i](r) - i]. 

^ i. rr rr' 


Previous applications of the CT-INT method to 
electron-phonon models only considered the case of opti¬ 
cal phonons. For the SSH model, the momentum depen¬ 
dence of the phonon dispersion leads to a pronounced 
minus sign problem whose origin can be traced back to 
the fact that r(r, r) (and hence the configuration weight) 
can take on negative values for r > 0 in the case of dis¬ 
persive phonons. The sign problem can be improved to 
some extent (at the expense of a higher expansion order) 
by generalized offsets a{r, s) (depending on Ising spins 
s = ±1, cf. Ref. [38]) in the bond operators B-(t) for 
vertices with negative weight. 


B. Observables 

We calculated the real-space correlation functions 


Ci,(r) = {(B^-h)(B^-h)) 

(bond) 

(12) 

Cp{r) = ((n^ -n)(no -n)) 

(charge) 

(13) 

CM = Mo) 

(spin) 

(14) 

CM = (AjAo) 

(pairing) 

(15) 


Here, n = {h-), b = {B-), and we consider the trans¬ 
verse spin correlations {{s^Sq) is identical within error 
bars to as required by symmetry). The pair¬ 
ing operator is defined as for spinful 

fermions (singlet pairing), and as A^ = for spin¬ 

less fermions (extended singlet pairing). 

In combination with the maximum entropy method 
[40], we calculated the single-particle spectral function 

A{k,o:) = I ^1 (*| 6, \j) , 

(16) 

with for spinless fermions, and for 

spinful fermions. Additionally, we calculated 

= ;^ E I - Cj), (17) 

with Og either a bond, charge, or spin operator. Here, 
|i) is an eigenstate with energy E-, and A^-- = E- — E-. 


IV. RESULTS 


The implementation of such interactions within the CT- 
INT method, including the use of auxiliary Ising fields to 
avoid a sign problem, has been discussed in Refs. [38, 39]. 
The simultaneous sampling of the interactions S-^ and S 2 
is done by randomly choosing a vertex type. 

The numerical effort scales as {a^L)^ with a ^ A for 
U = V = 0, and a - U for U > V, A. 


A. SSH model vs. optical SSH model 

Several different SSH-type models have been studied in 
the literature. The phase diagram for a half-filled band 
has been investigated numerically in terms of the origi¬ 
nal Hamiltonian (2) in Ref. [6] (with a cutoff to avoid an 
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FIG. 1. (Color online) Comparison between the optical SSH 
model (filled symbols) and the SSH model (open symbols) for 
(a) real-space bond correlations, (b) real-space spin correla¬ 
tions, (c) the single-particle Green function G{k = 0,t) and 
the dynamic bond correlation function AS'b(^ = 7r,T). Here, 
I3t = L = 22, cjJq = uj^ = O.lt, = y = 0, and A = 0.2. 
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FIG. 2. (Color online) (a) Bond, (b) charge, and (c) pairing 
correlations of the spinless optical SSH model as a function 
of ^ = Lsin(7rr/L) (see text). Here, /3t = L = 50, and 
ujQ — 0.5t. The solid lines correspond to 1/r^, the dashed line 
corresponds to 1/r. Here and in Figs. 6 and 11, the correlation 
functions are plotted only for odd distances r. 


discussed in Sec. V, we avoided the sign problem by sim¬ 
ulating the corresponding sign-free optical SSH models. 


unphysical sign change of the hopping term for large dis¬ 
tortions), in Ref. [15] by solving the original SSH model, 
in Ref. [25] by considering the optical SSH model (3), 
and in Ref. [31] by simulating the original interaction 
term in combination with a general phonon dispersion 
and extrapolating to acoustic and optical phonons [31]. 

Here, we present a comparison of the SSH model (2) 
and the optical SSH model (3), taking = O.lt to 

alleviate the sign problem; a larger value cJq = = 0.5t 

was used for most other results. Figure 1(a) shows the 
bond correlation function C^(r), whereas Fig. 1(b) shows 
the corresponding spin correlations. Remarkably, the re¬ 
sults for the models (2) and (3) are identical within sta¬ 
tistical errors. This quantitative agreement also extends 
to dynamic correlation functions, as demonstrated for 
the single-particle Green function at momentum /c = 0 
and the dynamic bond correlation function at g' = tt in 
Fig. 1(c). In addition to the results shown in Fig. 1, 
we compared bond, charge, spin, and pairing correlation 
functions, as well as the momentum distribution func¬ 
tion, for several other parameter sets, including nonzero 
U and y, higher temperatures, stronger coupling, and 
spinless models. All results were identical within error 
bars to those for the original SSH coupling and acoustic 
phonons. Based on this agreement, whose origin will be 


B. Spinless fermions 


Given the theoretical prediction of an extended metal¬ 
lic phase without long-range bond order in the spin¬ 
less SSH model [6], we show in Fig. 2 the real-space 
bond, charge, and pairing correlation functions defined 
in Sec. III. To eliminate boundary effects, we use the 
conformal distance ^ = L sin (Trr/L) [41]. This allows 
us to directly compare to the bosonization results for a 
spinless Luttinger liquid [42, 43] 




■ cos(2/cpr), 


CM) = 


K„ 


27r2p2 


+ 


■ cos(2/cpr), 


(18) 


For A = 0, all correlators decay with a power law 1/r^. 
According to our results, bond correlations become dom¬ 
inant for A > 0, suggesting < 1 and hence a repul¬ 
sive Luttinger liquid, in accordance with Ref. [17]. At 
the same time, charge and pairing correlations are sup¬ 
pressed. While the decay of pairing correlations can be 
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FIG. 4. (Color online) Single-particle spectral function 
A{k,uj) of the spinless optical SSH model. The dashed line 
indicates the Fermi level. Here, [31 = L = 30, ujq = 0.5t. 
Color scheme for this and other density plots from Ref. [47]. 


FIG. 3. (Color online) (a) Finite-size scaling of the bond 
correlations at the largest distance r = T/2, using a second- 
order polynomial extrapolation, (b) Finite-size scaling of the 
bond susceptibility for the spinless optical SSH model. Here, 
= L, and ujq — 0.5t. The key in (b) applies to both panels. 


explained by the increase of the exponent 2A:p 

charge correlations formally have the same exponent as 
2A;p bond correlations, see Eq. (18). The absence of en¬ 
hanced charge correlations in Fig. 2(b) can be attributed 
to a suppression of the amplitude, ^ 0 [44], causing 
the first term (1/r^) to dominate in the metallic phase 
[the exponent is close to 2 in Fig. 2(b)]. Pairing corre¬ 
lations [Fig. 2(c)] are always subdominant in the metal¬ 
lic phase. In the insulating phase, charge and pairing 
correlations are strongly suppressed—consistent with an 
exponential decay—while bond correlations become long- 
ranged. 

To determine the critical value for the Peierls tran¬ 
sition, we show in Fig. 3(a) a finite-size scaling of the 
bond correlations at the largest distance r = L/2. The 
extrapolation using system sizes up to L = 50 suggests 
= 0.9(1). The onset of long-range order can also be 
tracked by calculating the bond susceptibility [25] 


(19) 


As a function of L (and for /3 = 1/), Xp(7r)/I/ is expected 
to vanish in the metallic phase, and to diverge ^ L in 
the dimerized Peierls phase. For spinless fermions, Xh 
has a more favorable scaling with system size than the 
bond correlations shown in Fig. 3(a), and does not involve 
an unknown scaling function. Our numerical results in 
Fig. 3(b) confirm the expected behavior of the suscepti¬ 
bility, and give a critical value consistent with Fig. 3(a). 
We also carried out simulations for other phonon frequen¬ 
cies, and found that A^ increases with increasing cJq, in 
agreement with theoretical expectations [6]. 

Within a bosonization description of the Luttinger liq¬ 


uid to (charge-density-wave) insulator transition, = 
1/2 corresponds to the critical point of the Kosterlitz- 
Thouless phase transition. Equation (18) therefore sug¬ 
gests a 1/r decay at the critical point. Although we have 
to consider the possibility of logarithmic corrections in 
the vicinity of the transition, the correlation functions in 
Eig. 2(a) indeed decay slightly faster than 1/r (consistent 
with > 1/2) for A = 0.8, and slightly slower than 1/r 
for A = 1, consistent with A^ = 0.9(1). 

To further characterize the two phases, we calculated 
spectral functions. Eigure 4 shows the single-particle 
spectral function for (a) A = 0.6 (metallic phase) and (b) 
A = 2 (insulating phase). In Eig. 4(a) we see a gapless 
(within the available momentum resolution) band derived 
from the bare dispersion —2t cos k. The most prominent 
effects of the electron-phonon coupling are (i) a substan¬ 
tial broadening outside the coherent interval [—cJq,cJq], 
and (ii) an enhanced bandwidth [^ 6t in Eig. 4(a)] arising 
from the enhanced hopping due to the SSH interaction, 
see Eq. (3). (The bandwidth of the spectrum remains 

4t in Holstein-type models [45].) In the Peierls phase, 
Eig. 4(b), the spectrum has acquired a gap at the Eermi 
level, and exhibits signatures of backfolded shadow bands 
characteristic of systems with reduced periodicity [46]. 
While these features are qualitatively captured by mean- 
field theory [45], at and near /cp, we can identify two sep¬ 
arate energy scales in Eig. 4(b) previously discussed for 
the spinless Holstein model [45]. The lowest-energy exci¬ 
tation can be related to soliton excitations [45], whereas 
the high-energy excitations correspond to the renormal¬ 
ized cosine band structure with a dimerization gap. 

Eigure 5 shows the corresponding dynamic bond and 
density structure factors. Eor weak coupling, the latter 
still resemble the noninteracting particle-hole continuum. 
Additionally, the bond structure factor in Eigs. 5(a) and 
(b) reveals at low energies the renormalized phonon fre¬ 
quency, in particular a central (cj = 0) peak at g = tt 
in the Peierls phase in Eig. 5(b) related to long-range 
order (in the Holstein model, where the lattice is con- 
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FIG. 6. (Color online) (a) Bond and (b) spin correlations 
of the spinful optical SSH model as a function of ^ = 
Lsin(7rr/L). Here, /3t = L = 50, and ujq — 0.5t. The solid 
lines correspond to 1/r^ 


FIG. 5. (Color online) (a)-(b) Dynamic bond and (c)-(d) 
dynamic density structure factor of the spinless optical SSH 
model. Here, = L = 30, and uq — 0.5t. 


pled to the charge, the phonon dispersion appears in the 
dynamic charge structure factor [45, 48]). The density 
structure factor in the metallic phase, Fig. 5(c), features 
a gapless linear mode at small g, whereas a gap and a 
strong suppression of low-energy excitations are visible 
in the Peierls phase, see Fig. 5(d). 


C. Spinful fermions 

We begin by considering the real-space correlation 
functions shown in Fig. 6. Similar to the spinless case, 
bond correlations decay increasingly slowly upon switch¬ 
ing on the electron-phonon coupling. At the same time, 
we observe a suppression of charge, pairing and spin cor¬ 
relations [spin correlations are shown in Fig. 6(b), the 
others behave similarly]. Since we expect an ordered 
Peierls state for any A > 0, bond correlations should 
oscillate with a constant, nonzero amplitude at large dis¬ 
tances, and all other correlation functions should decay 
exponentially. While the finite-size extrapolation of the 
bond correlations at the largest distance in Fig. 7(a) gives 
a nonzero value already for A = 0.2, the expected expo¬ 
nential decay of spin correlations is only clearly visible 
(for the system size used) for A = 0.6 in Fig. 6(b). 

The Peierls state can be detected even at weak cou¬ 
pling. Because the bond susceptibility is not a reliable 
indicator to track the onset of long-range order in systems 
with a spin gap [49], we instead consider the finite-size 
estimates of the Luttinger liquid parameters and K^. 
The latter are defined in terms of the charge and spin 


structure factors as 

K^iL) = nS^iq,)/q, , 

K^{L) = nSM/Qi . ( 20 ) 

where q-^ = 27 t/L is the smallest nonzero wave vector for 
system size L. Finite-size extrapolation has been shown 
to give accurate results for the Luttinger parameters, for 
example, in the extended Hubbard model [50]. However, 
as argued in Refs. [25, 51], even very small spin gaps 
can be detected from values K^{L) < 1, as they imply 
K^{L oc) = 0. In contrast, in systems with gap¬ 
less spin excitations, numerical simulations usually give 
KW > 1 with a very slow (logarithmic) convergence to 
the value = 1 implied by the SU{2) spin symmetry. 
This diagnostic gives reliable results for the integrable at¬ 
tractive and repulsive Hubbard models [52]. For the SSH 
model, we find K^{L) < 1 and hence a spin gap even 
for A = 0.1 [Fig. 7(b)], which is consistent with a Peierls 
state for any A > 0. The same behavior is observed over 
the whole range of phonon frequencies from the adiabatic 
to the antiadiabatic regime, see Fig. 7(c). 

In general, a spin gap does not necessarily imply long- 
range Peierls order. For example in the Holstein model, 
an extended metallic (Luther-Emery [53]) phase with a 
spin gap but gapless charge excitations [54] exists. While 
theory suggests that the spin and charge gaps open si¬ 
multaneously in the SSH model, we can obtain explicit 
evidence from the degeneracy of K^{L) and K^{L) in 
Fig. 7(b). Within statistical errors, the correlation func¬ 
tions Sp{q^T) and S^{q^T) and hence also K^{L) and 
K^{L) are identical. This is a consequence of the spin- 
charge symmetry discussed in Sec. II, and has impor¬ 
tant consequences: if K^{L) < 1 indicates a spin gap, 
™plies also a charge gap and hence 
insulating behavior. Figure 7(b) therefore suggests that 
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FIG. 7. (Color online) (a) Finite-size scaling of the bond 
correlations at the largest distance r = L/2 for the spinful 
optical SSH model and different values of A. Here, = L, 
and ujo = 0.5t. The inset shows a closeup of the extrapolated 
values, (b) Finite-size estimates of the Luttinger liquid pa¬ 
rameters [we find Kp(L) = Ka-{L) within error bars] for the 
same parameters as in (a), (c) Kp(L) and K(j{L) for different 
phonon frequencies. Here, fit — and A = 0.1. 
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FIG. 8. (Color online) Single-particle spectral function 
A{k,uj) of the spinful optical SSH model. The dashed line 
indicates the Fermi level. Here, fit = L = 30, and c^o = 0.5t. 
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FIG. 9. (Color online) (a) Dynamic bond and (b) dynamic 
spin structure factor of the spinful optical SSH model. Here, 
fit = L = 30, ujo = 0.5t, and A = 1. 


the spinful optical SSH model is a Peierls insulator even 
for very small A, in agreement with the theoretical pre¬ 
diction of A^ = 0. Given the quantitative agreement be¬ 
tween the optical SSH model and the SSH model demon¬ 
strated in Fig. 1, our findings carry over to Eq. (2). 

Having established the insulating behavior for any 
A > 0, we now consider excitation spectra. Figure 8 
shows the single-particle spectral function for A = 0.6 
and A = 1. The gap opens exponentially starting from 
A = 0. Nevertheless, we see a small gap and backfolded 
shadow bands in Fig. 8(a), in contrast to the results for 
the metallic phase of the spinless model in Fig. 4(a) for 
the same value of the coupling (A = 0.6). For a stronger 
coupling A = 1, Fig. 8(b) reveals a substantial gap, and 
the spectrum closely resembles that of the spinless SSH 
model in the Peierls phase, see Fig. 4(b). 

Figure 9 shows the dynamic bond and spin struc¬ 
ture factors for A = 1. As expected for the dimerized 
phase, Fig. 9(a) features a central peak with large spec¬ 
tral weight 8it q = TT. Both charge and spin excitations 
(identical within error bars) are gapped at small q. 


D. Spinful fermions with Coulomb interaction 

We now turn to the spinful SSH-t/E model with U = 
AV = 2.5t, considered to be appropriate for polyacetylene 
[7], and a phonon frequency cJq = 0.5t. In principle, our 
method can also be applied to models with long-range 
interactions [55]. Because Eqs. (2) and (3) give identical 
results even with the additional electron-electron inter¬ 
action, we simulated the optical SSH model to avoid the 
minus-sign problem mentioned in Sec. III. 

Erom Refs. [25, 31, 56], it is known that the phase dia¬ 
gram consists of Mott and Peierls insulating regions, ex¬ 
cept for the adiabatic limit where the Peierls state seems 
to remain stable for any A > 0. Away from this limit, 
the system is a Mott insulator for A < A^, with gapless 
spin excitations but a gap for single-particle and charge 
excitations. Eor A = 0 (extended Hubbard model), the 
ground state is a Mott insulator for U ^ 2V [44, 51, 57], 
as relevant for polyacetylene. The electron-phonon cou¬ 
pling competes with the dominant antiferromagnetic spin 
correlations, and eventually gives rise to a Peierls phase 
with a spin gap and long-range bond order [25, 31]. 

To estimate A^ for the Mott-Peierls transition for our 
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FIG. 10. (Color online) Finite-size scaling of Ka(L) for the 
optical SSH-f/y model. Here, = L, uo = 0.5t, U = 2.5t, 
and V = 0.625t. 


choice of parameters, we carried out a finite-size scaling 
of K^{L); the additional Coulomb interaction restricts 
us to system sizes L < 30. The results in Fig. 10 for 
A = 0.1 are consistent with the absence of a spin gap and 
= 1, and therefore with a Mott insulating state. In 
contrast, for A > 0.2, we have K^{L) < 1 and a finite-size 
dependence consistent with = 0 and a spin-gapped 
Peierls phase. Moreover, we find K^{L) < 1 for all values 
of A considered, and a very weak dependence of K^{L) 
on A in the range 0 < A < 0.4 (not shown). 

Figure 11 shows the real-space correlation functions 
for bond, spin, charge and pairing operators in the Mott 
(A = 0.1) and the Peierls phase (A = 0.4), respectively. 
Because the system has a charge gap in both phases, the 
charge and pairing correlations (see Ref. [42] for the ex¬ 
plicit Luttinger liquid results) in general decay exponen¬ 
tially; the former because = 0 and ^ 0, similar 
to the spinless case, the latter because the exponents di¬ 
verge for = 0. Figures 11(c) and 11(d) reveal a strong 
suppression compatible with an exponential decay. For 
nonzero U and R, the SSH-t/R model is no longer invari¬ 
ant under the modified particle-hole transformation, and 
the degeneracy of spin and charge correlations is lifted. 

The bond and spin correlations are more interesting. 
In a Luttinger liquid, their 2k^ components decay as [58] 

cf'(r) = COS(2t,r), (21) 

where we neglected potential logarithmic corrections due 
to marginal operators. In the Mott phase with = 0 
and = 1 (corresponding to a Luttinger liquid in the 
spin sector, or COSl in the notation of Ref. [59]), Eq. (21) 
predicts a 1/r decay, in good agreement with our results 
in Fig. 11(a) and 11(b). While we cannot exclude log¬ 
arithmic corrections on the accessible system sizes, our 
results are consistent with a larger amplitude for spin 
correlations than for bond correlations, > ^b* 

Peierls phase = 0, COSO), bond correlations 

become long-ranged and spin correlations decay exponen¬ 
tially as a result of the nonzero spin gap, see Fig. 11. 

Although finite values of U, V give rise to a Mott 
phase for small enough A, it has been observed that 2/cp 
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FIG. 11. (Golor online) (a) Bond, (b) spin, (c) charge, and (d) 
pairing correlations of the optical SSH-f/E model as a func¬ 
tion of ^ = Lsin(7rr/L). The solid (dashed) line illustrates a 
1/r^ (a 1/r) decay. Here, /3t = L = 30, c^o = 0.5t, U = 2.5t, 
and V = 0.625t. 


bond correlations are actually enhanced by Coulomb in¬ 
teractions [23, 24, 34, 60]. Our results in Fig. 11(a) con¬ 
firm such an enhancement. Interestingly, previous work 
showed that the power-law exponent of 2k^ bond corre¬ 
lations is independent of U and V [24], namely equal to 
one in the Mott phase and zero in the Peierls phase; the 
data in Fig. 11(a) are compatible with this result. The 
enhancement of bond correlations is therefore related to a 
larger amplitude A^. While A^ grows continuously with 
increasing U, A^ has been found to have a maximum near 
l/-4t [23, 24, 34, 56]. 

Figure 12 shows the single-particle spectral function 
for A = 0.1 [Mott phase, Fig. 12(a)] and A = 0.4 [Peierls 
phase, Fig. 12(b)]. In both cases, we find a gap at the 
Fermi level that increases with increasing A. Moreover, 
the SSH electron-phonon coupling again gives rise to an 
enhanced bandwidth of the cosine-like band. For in¬ 
termediate values near the Mott-Peierls transition (not 
shown), we find no evidence for a closing of the gap, con¬ 
sistent with an insulator-insulator transition without an 
intermediate metallic state. 
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FIG. 12. (Color online) Single-particle spectral function 
A{k,uj) of the optical SSH-f/l/ model. The dashed line indi¬ 
cates the Fermi level. Here, — L — 30, c^o = 0.5t, U = 2.5t, 
and V = 0.625t. 


In Fig. 13, we present the two-particle excitation spec¬ 
tra in the Mott phase (A = 0.1) and the Peierls phase 
(A = 0.4). Whereas only a partial softening of the phonon 
mode 8it q = TT is visible in the structure factor in the 
Mott phase [Fig. 13(a)], we observe a central peak at 
g = TT in the Peierls phase [Fig. 13(b)]. The dynamic 
density structure factor reveals a gap for long-wavelength 
excitations in both cases, see Figs. 13(c) and 13(d). Fi¬ 
nally, the dynamic spin structure factor in the Mott phase 
[Fig. 13(e)] has a gapless linear mode at small g, and en¬ 
hanced spectral weight at g = tt = 2k^ reflecting the 1/r 
antiferromagnetic spin correlations. In the Peierls phase, 
Fig. 13(f), the weight at g = tt is suppressed, but the 
spin gap at g = 0 is much smaller than in Fig. 9(b). 


V. DISCUSSION 



0 7r/2 TT 0 7r/2 tt 



0 7 t /2 tt 0 7 t /2 tt 



0 7 t /2 tt 0 tt/2 tt 


Q q 

FIG. 13. (Golor online) (a)-(b) Dynamic bond, (c)-(d) dy¬ 
namic charge, and (e)-(f) dynamic spin structure factor of 
the optical SSH-UP model. Here, /3t = L = 30, c^o = 0.5t, 
U = 2.5t, and V — 0.625t. 


A. Agreement of SSH and optical SSH models 

Figure 1 reveals a remarkable quantitative agreement 
between results obtained for the SSH model (2) and 
the optical SSH model (3). We have found this agree¬ 
ment to hold for all parameters and observables con¬ 
sidered. While theoretical arguments suggest that only 
the phonon frequency at g = 2/cp is important for the 
low-energy physics of the SSH model [6, 16, 56], devi¬ 
ations may be expected at higher energies or smaller 
length scales, as typically probed in numerical simula¬ 
tions. However, for the parameter sets considered, the 
results are identical within error bars. This equivalence 
of the models (2) and (3) is consistent with the excellent 
agreement of the phase boundary for the optical SSH-I/U 
model from quantum Monte Carlo simulations [25] and 
that obtained with the functional renormalization group 
[18] for the SSH-I/U model with acoustic phonons. On 
the other hand, a significant shift of the critical value 
for the Mott-Peierls transition was reported in Ref. [31], 
although the values become almost identical for small U 
and V; however, these results were obtained for a model 


with the original SSH coupling term. Therefore, the op¬ 
tical limit does not correspond to the optical SSH model 
considered here and in Ref. [25]. 

To understand the quantitative agreement between the 
SSH model and the optical SSH model, it is useful to con¬ 
sider the relation of the two models. The variables Q- of 
the SSH model (2) describe the longitudinal displacement 
of lattice site i from its original position. The hopping 
integral between sites i and i +1 depends on the length of 
the bond between these two sites, and hence on ~Qi' 

A change of also changes the bond between sites i and 
i — 1 because the lattice displacements at neighboring 
sites are coupled; this coupling gives rise to an acoustic 
phonon mode. In contrast, in the optical SSH model (3), 
the variables describe the bond length itself, that is, 
they correspond to the difference ~ Qi Eq. (2). 
Each of the L bonds is modeled as a harmonic oscillator 
with displacement Q- and momentum P-. Neighboring 
bonds are considered to be independent, as described by 
Einstein phonons. While the impact of this approxima¬ 
tion is a priori unclear, our numerical results suggest that 
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both models give identical results. 

We attribute the quantitative agreement to the fact 
that at half-filling (2/cp = tt), the 2/cp bond correlations— 
generic for Luttinger liquids and Mott insulators [with 
a power-law decay as described by Eqs. (18) and (21)] 
as well as for Peierls insulators with long-range bond 
order—are equivalent to alternating displacements of 
neighboring bonds. Because of the inherent preference 
for 2/cp (alternating) bond correlations, it is not neces¬ 
sary to encode the anticorrelation between the lengths 
of neighboring bonds into the Hamiltonian. While the 
SSH and optical SSH models are not connected by a 
canonical transformation, a change to bond variables 

~ (Qi +1 ~ Qi)!^ conjugated bond momenta p- 
in Eq. (2) gives the optical SSH model (3) up to an ad¬ 
ditional term of the form responsible for the 

phonon dispersion. Neglecting this term amounts to re¬ 
placing = co’^sin(g/2) by Because sin(g/2) 

has a maximum at g = tt, the first-order correction ^ q 
vanishes for half-filling at the dominant value q = 2k^ 
that determines the low-energy physics. Erom both argu¬ 
ments provided, we expect the models (2) and (3) to give 
different results away from half-filling where 2k^ ^ tt. 


B. Lattice fluctuations and low-energy theory 

The question of the impact of quantum lattice fluctua¬ 
tions on the dimerized Peierls state has been raised soon 
after the original work of SSH [6, 11, 12]. However, while 
(thermal or quantum) fluctuations are in principle ex¬ 
pected to destroy long-range order for sufficiently weak 
electron-phonon interactions, the Peierls state remains 
stable at T = 0 in the spinful SSH model (although the 
dimerization is reduced), as confirmed by renormaliza¬ 
tion group methods [17, 18, 56, 61] and numerical simu¬ 
lations [6, 15, 25, 31]. In this respect, the SSH model is 
hence quite different from the Holstein model. In the lat¬ 
ter, a metallic phase emerges from lattice fluctuations in 
both the spinless and the spinful case [19-21]. The spin¬ 
ful SSH model hence constitutes a case where Peierls’ 
theorem is valid beyond the adiabatic limit. 

Eor the SSH model, lattice fluctuations destroy long- 
range order if additional electron-electron interactions 
are taken into account [25, 31]. In this case, the criti¬ 
cal value depends strongly on the phonon frequency 
[25]. However, a metallic state does not seem to exist. 
Instead, the phase diagram in the (A,co’q) plane for fixed 
large U reflects the physics of the Heisenberg spin-phonon 
model, which is always dimerized in the adiabatic limit, 
but has a gapless phase with critical spin correlations for 
cJq > 0 and A < A^ [25, 28, 62, 63]. Eor large t/, this 
phase corresponds to the Mott phase of the SSH model. 

The absence of a metallic phase in the spinful SSH 
model was also predicted from two-cutoff renormaliza¬ 
tion group calculations [17, 56, 61]. Because the forward 
scattering matrix element ^2 is exactly zero, any repul¬ 
sive umklapp scattering is relevant. While this result 


agrees with numerical results and the functional renor¬ 
malization group, applications of the same method to 
the Holstein model [17, 56, 61] fail to describe the ex¬ 
tended metallic phase at weak coupling. The two-cutoff 
renormalization group is hence not reliable in general. 
The likely reason why it gives correct results for the SSH 
model is that the physics of the latter at finite uJq is es¬ 
sentially the same as for cJq = oo [6], so that retardation 
effects do not play a significant role. 

While a purely fermionic theory is in general not suf¬ 
ficient to capture the physics of systems with a coupling 
to quantum phonons, a comparison with bosonization 
and renormalization group results for fermionic models 
provides valuable insights. Although the microscopic 
problem is rather different, the metal-insulator transi¬ 
tion of the spinless SSH model can be related to a spin¬ 
less Luttinger liquid that undergoes a Mott transition to 
a gapped 2k^ charge-density-wave phase. In fact, in the 
limit cJq = oo, the spinless SSH model maps to a model of 
spinless fermions with hopping t and nearest-neighbor re¬ 
pulsion V = X, which exhibits a quantum phase transition 
from a Luttinger liquid to a charge-density-wave insula¬ 
tor ai V = 2t (corresponding to = 1/2) [43]. Similarly, 
numerical results for the spinless Holstein model—for 
which a strong-coupling picture also leads to a model of 
spinless fermions with repulsion v —also support a metal- 
insulator transition at = 1/2 [64]. 

The physics of the spinful SSH and SSH-l/E models 
may be related to bosonization results for the extended 
Hubbard model with onsite interaction u and nearest- 
neighbor interaction v [44]. Here, u and v correspond 
to effective interactions that emerge from both electron- 
electron and electron-phonon interactions. Eor the case 
U = V = 0, the existence of spin and charge gaps for any 
A > 0 may be related to the fact that for = 0, umklapp 
and backscattering are relevant for any > 0, giving rise 
to a fully gapped phase with long-range 2kp charge cor¬ 
relations (rather than a bond-order wave as in the SSH 
model; for the extended Hubbard model, 2kp bond corre¬ 
lations have the same exponent as 2kp charge correlations 
but a vanishing prefactor [44]). On the other hand, the 
Mott-Peierls transition observed for nonzero Coulomb in¬ 
teractions 1/ and V (with U V) may be understood 
in terms of the extended Hubbard model with > 0 and 
V > 0. Eor the latter, umklapp scattering is relevant for 
any u > 0 and gaps out the charge sector [44], while 
a transition to a spin-gapped 2/cp charge-density-wave 
phase takes place near u = 2v [44]. The Mott phase 
with gapless spin excitations that exists for u ^ 2v is 
identical to the Mott phase observed for A < A^ in the 
SSH-l/E model, whereas the fully gapped charge-density- 
wave phase may be related to the dimerized Peierls phase. 
(Eor the extended Hubbard model at weak to intermedi¬ 
ate couplings u and v, two independent transitions for 
the spin and charge excitations as well as an intermedi¬ 
ate bond-ordered phase exist [51, 57, 65-68].) 

Although a metallic phase is absent in the SSH model, 
we note that an extension of an argument by Voit for 
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charge-density-wave insulators [69] implies that a metal¬ 
lic phase with dominant bond correlations can only be 
realized in Luther-Emery liquids with a gap for spin ex¬ 
citations, because according to Eq. (21) spin and bond 
correlations decay with the same power-law exponent in 
a Luttinger liquid with = 1. 

Although the bosonization results for the extended 
Hubbard model capture several key aspects of the SSH- 
UV model, the continuum field theory does not account 
for the spontaneous breaking of a discrete Ising Z 2 sym¬ 
metry at the Peierls transition. Instead, it provides a 
theory of U{1) charge and spin phase fields that can be 
pinned when umklapp and/or backscattering terms be¬ 
come relevant, leading to Kosterlitz-Thouless transitions 
and massive excitations but no long-range order. Ac¬ 
cordingly, the Mott phase has a charge gap and critical 
spin correlations, but does not break any symmetries. In 
contrast, the Peierls phase is characterized by charge and 
spin gaps and a broken Ising symmetry reflected in the 
central peak visible in Eig. 13(b). The existence of a 
U{1) X symmetry in the continuum limit has recently 
been illustrated for the spinless Holstein model [48]. 

VI. SUMMARY 

We used the CT-INT quantum Monte Carlo method 
to study SSH models with quantum phonons. Eirst, we 
demonstrated that the original SSH model with acoustic 
phonons gives results that agree within statistical errors 
with those for the optical SSH model, a property that 


can be related to the inherent tendency toward 2k^ bond 
correlations. Next, we studied the Peierls metal-insulator 
transition of the spinless SSH model in terms of real-space 
correlation functions and excitation spectra. The results 
are consistent with a Kosterlitz-Thouless quantum phase 
transition at a critical Luttinger parameter = 1/2. 
The spectra reveal the expected Peierls gap, soliton sig¬ 
natures, and a central peak related to long-range order. 
In the spinful case, we confirmed the absence of a metallic 
phase by detecting a spin gap even at weak coupling, and 
by demonstrating a symmetry-related degeneracy of spin 
and charge correlations. The spectral functions reveal a 
gap for single-particle, charge and spin excitations, and a 
central peak. Einally, we studied the spinful SSH model 
with additional electron-electron interactions appropri¬ 
ate for polyacetylene. Our results are consistent with a 
phase transition from a Mott state with critical spin and 
bond correlations to a Peierls state with long-range bond 
order and exponential decay of all other correlation func¬ 
tions. Electron-electron repulsion was found to enhance 
bond correlations. In contrast to the Holstein-Hubbard 
model, we found no indications for metallic behavior. 
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